Entanglement of Collectively Interacting Harmonic Chains: An effective 

Two-Dimensional System 
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We study the ground-state entanglement of one-dimensional harmonic chains that are coupled 
to each other by a collective interaction as realized e.g. in an anisotropic ion crystal. Due to the 
collective type of coupling, where each chain interacts with every other one in the same way, the total 
system shows critical behavior in the direction orthogonal to the chains while the isolated harmonic 
chains can be gapped and non-critical. We derive lower and most importantly upper bounds for the 
entanglement, quantified by the von Neumann entropy, between a compact block of oscillators and 
its environment. For sufficiently large size of the subsystems the bounds coincide and show that the 
area law for entanglement is violated by a logarithmic correction. 



Presently there is a growing interest in the interrela- 
tion between entanglement and ground-state properties 
of many-body lattice models. For a number of spin sys- 
tems [l[ a strict correspondence between the absence of 
criticality, the presence of an energy gap, and an area 
law for the entanglement was established. The latter 
states that the entanglement of a compact sub-set of lat- 
tice sites with the rest of the system, measured by the 
von Neumann entropy, scales with the surface area of the 
sub-set. For critical spin systems it was shown that an 
additional logarithmic correction to the area law emerges. 
A similar relation between criticality and entanglement 
was suggested for harmonic lattice models 0, HJ . In [1, [f| 
an area law was established for harmonic lattice models 
in arbitrary dimensions with nearest-neighbor coupling 
which have a gaped spectrum. For finite-range couplings 
in one dimension a one-to-one correspondence between 
the validity of the area law and non-criticality was estab- 
lished in [fj , and logarithmic corrections were derived for 
critical systems. 

Although the relation between criticality and entropy- 
area law seems rather universal, there are a number of 
examples where this relation does not hold [H, 0] ■ Until 
now there is no general understanding of the conditions 
for the validity of an entropy area law in particular in 
higher dimensions 0, 0, Q • In the present paper we 
discuss a specific gapless oscillator model with dimension 
larger than one, for which an exact asymptotic expres- 
sion for the entropy can be obtained. Due to the col- 
lective nature of the interactions in one spatial direction 
the system is critical and thus a violation of the area law 
is expected. We here derive a lower and, most impor- 
tantly, a tight upper bound for the entropy and in this 
way obtain an exact form of the correction term to the 
area law. 

Let us consider a set of parallel harmonic chains (see 
Fig[T]) each containing n x oscillators, with n x — > oo in the 
thermodynamic limit. We will refer to the direction par- 
allel to the chains as x-axis, and to the orthogonal direc- 
tion as j/-axis. The number of parallel chains is denoted 
as n y , again with n y — > oo in the thermodynamic limit. 
The oscillators are described by the canonical variables 



(<Zi>Pi)j where i = 1,2, ...N (N — n x n y ) is a collective 
index that labels the oscillator. We adopt the following 
notation: i = 1, n x correspond to the oscillators in the 
first chain with growing x coordinate, i = n x + 1, 2n x 
corresponds to oscillators in the second chain and so on. 
We consider a quadratic Hamiltonian of the form 



1 N 1 N 



(1) 



with a coupling matrix V. We are interested only in 
a translationaly invariant coupling, i.e. we assume that 
the matrix elements of V depend only on the difference 
of the x coordinates and the difference of the y coordi- 
nates. Hence V is a block Toeplitz matrix. For oscillator 
systems with a quadratic coupling of the form of eq. ([T]) 
the ground state 



*o (q) ~ exp 



I <q| V 1 ' 2 |qj 



(2) 



and all its properties, as e.g. the correlation length in po- 
sition or momentum space, are determined by the square 
root of V, where q = (qi , qi, . . . , qjq) is the vector of posi- 
tion variables. The ground state can easily be determined 
if V is the square of another matrix, which we assume to 
be again a Toeplitz matrix, 



V = Z 2 /n y . 



(3) 



The factor l/n y is choosen such that the matrix elements 
of V remain finite in the limit N — > oo. Assuming Z to 
be a Toeplitz matrix guarantees that the coupling V is 
a Toeplitz matrix as well. We furthermore consider Z to 
be of the block-matrix form 



(4) 



The elements of Z are n x x n x matrices and characterize 
according to eq.© correlations. The diagonal elements 
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of Z describe correlations within one chain, i.e. in x 
direction, the off-diagonal elements describe correlations 
between the chains. A and Q are both assumed to be 
Toeplitz matrices of finite range, i.e. their matrix ele- 
ments Kk and Qk, where A& = Aj_u_j| = (i|A|j), vanish 
exactly for k > R. The finite range of A and Q ensures 
that the interaction V is of finite range within the chains, 
while the form of Z implies that V is constant orthogo- 
nal to the chains. We assume furthermore that A, Q and 
A — Q are positive definite matrices. A simple calculation 
shows that the ground state of V is degenerate and in the 
thermodynamic limit n x , n y — > oo has only one non-zero 
eigenvalue. This means that the total Hamiltonian, Eq. 
(HJ), is gapless. It should be noted however that the col- 
lective nature of the interactions is not sufficient for a 
gapless spectrum of the Hamiltonian. 

Since all off-diagonal elements of Z are identical, cor- 
relations between oscillators do not depend on their dis- 
tance in y direction and the total system is critical irre- 
spective of the correlation properties within the chains. 
Thus one expects that the entropy area law is broken. 
In fact one can easily find a lower bound to the en- 
tropy by the following simple argument: Let us consider 
a partition of the set of N oscillators into a compact 
sub-system I with No = l x l y and a sub-system II with 
N — No oscillators (see FigJTJ. If we now consider har- 
monic chains in y— rather than in x— direction, the "y"- 
chains couple to each other with finite-range interaction 
A (see Fig[T] b). We thus have reason to assume that 
S > I x Sq, where Sq is the entropy of a single "y" -chain. 
Since the coupling within the chain is now collective (Q), 
the "y" -chain itself is critical and its entropy scales as 
So ~ Inly. Thus S > l x In l y which in the thermody- 
namic limit {l x ,l y } - * oo is much larger than surface 
area 2(l x + l y ). While it is easy to see that the area law 
is broken, it is non-trivial to find an upper bound to the 
entropy and the exact form of the correction term. This 
will be done in the following. 



matrices V 1 ! 2 and V can be decomposed as 



V 1/2 = [(A - Q) <g> l y + n v Q <g> V ny ,n y ] /y/ni, (5) 



and 



v- 1 ' 2 = | (A - Qy 1 ®ly + 
■[(A - Q + nyQ)- 1 - (A - Q)- 1 ] <g) V„ y , ny ], 



(6) 



where l y is the unity matrix of size n y x n y and V nm — 
\Pnm)(Pnm\ is the projector onto the (in general non- 
normalized) vector 



Following Refs. @, HI, 0, Q , the von-Neumann entropy 
or the entropy of entanglement of the two compact parts 
I and II can be calculated from a decomposition of V 1 ^ 2 
into the two subsystems. To this end we express V 1 / 2 and 
y-i/ 2 i n a block form according to the two sub-systems 
by proper reordering of rows and columns 



y-V2 = 



A B 
B T C 



V 1 ' 2 = 



D E 
E T F 



(7) 



Here A and D are Ao x No matrices describing corre- 
lations within sub-system I, C and F are (N — No) x 
(N — No) matrices describing correlations within sub- 
system II, and the matrices B and E describe the cor- 
relations between them. The entropy of entanglement 
is then given by the eigenvalues fa > 1 of the matrix 
product A ■ D 



N 



fix) 



i=l 
X + l 



In- 



(8) 



(9) 



(a) 



• 








• 




m 


• 


e 




Q 


• 


Q 


• 




d 






• 




• 










• 




• 










• 




• 


• 








• 




• 


• 








• 




• 







• ) 




• } 


C f 




• ) 


c * 




• ) 


c • 




• ) 


o 




• ") 


r • • • 





(b) 



Despite the simplicity of its form, expression ([5]) cannot 
be explicitly evaluated in general. Due to the special 
interaction matrix the eigenvalues can however be evalu- 
ated in the thermodynamic limit: 

From the spectral decomposition of V 1 / 2 , eq.©, one 
easily finds that the subsystem matrices read 



A = [A Q <x> h y + (Ai - Ac) <8> V ny ,i y ] 



(10) 



D=\D Q <E> h v + n v D 1 ® V n 



FIG. 1: (a) Collectively interacting strings of harmonic os- 
cillators with finite-range intra-chain coupling A and collec- 
tive inter-chain coupling Q. The grey area indicates the 
sub-system I of oscillators, (b) Alternative view: interacting 
strings with collective intra-chain coupling Q and finite-range 
inter-chain coupling A. 

Using the spectral representation of V, the correlation 



where Aq,Ai and Dq,Di are l x x l x principal subma- 



trices of (A - Q) _1 , (A - Q + n y Q)~ 
respectively. For large n y one has 



and (A - Q), Q 



A ■ D » (Ao ■ Do) ® h y + n y (A ■ Di) ® V ny ,i y ■ (11) 

Here we have used that V 2 , — l v /n v V n „ i„ which scales 
as l/n y for fixed l y and is thus negligible in the ther- 
modynamic limit. Furthermore V n i has one nonzero 



3 



eigenvalue l y /n y , which vanishes in the thermodynamic 
limit (l y fixed and n y — ► oo), and (L — 1) zero eigenval- 
ues. Thus the l x l y eigenvalues of A-D can be decomposed 
into two sets. The first set consists of the l x eigenvalues 
of Aq ■ Dq each of which occurs (l y — 1) times: 

Mi, ' ' ' = a i ( A o ■ D ) , 

Uly,--' !^2(i v -i) = a 2 (A -D ), (12) 



M({ e -i)(J y -i)+i) ' ' ' >Mi !D (i s ,-i) - (A) • A)) • 

Here and in the following (X) denotes the fcth eigen- 
values of the matrix X. The total number of these eigen- 
values is l x (ly — 1)- The second set consists of the l x 
eigenvalues of (A Q ■ D + l y (A ■ D±)) 



li k = u k (A ■ D + l y (A • £>i)) , 
for k = l x (l v - 1) + 1, 



(13) 



■7 lxly- 



Expression ([13"]) for the second set of eigenvalues can be 
simplified using Lidskii's theorem [l(| which states: Let 
X and Y be M -dimensional Hcrmitian matrices. More- 
over let oifc (X) ,aik (Y) and a% (X — Y) , k = 1, ...,M be 
the eigenvalues of X, Y and X — Y respectively in as- 
cending order {ct\ (X) < a 2 (X) < ... < ctu (X)}. Then 
there exist numbers Wkj > 0, (k,j = 1, M), such that 
Efe w kj = Ej ^fej = 1 and 



a fe (X) = a k (Y) + J2 «OT (X - Y) . 



(14) 



Equation ([14]) implies that for sufficiently large Z w the 
eigenvalues of the matrix Aq ■ Do + l y (Aq ■ D%) are 



a k (A ■ D + ly (A ■ D-j)) w l y }^ w kj aj (A ■ Di) 

j'=i 



An upper bound to the entropy can be found by evalu- 
ating the sum over the eigenvalues ([12")) and (TT3]) in eq. ([5]) 
separately 



S = S1+S2 



(16) 



E 

3=1 



Fl 



l x ly 

E 

3'='aj(il/-l) 



Taking into account eq. ([T2]) one recognizes that 5i is 
apart from a prefactor (l y — 1) formally equivalent to 
the von-Neumann entropy of a linear oscillator chain of 
length l x with interaction V = (A — Q) 2 



entropy of such a linear chain saturates in the thermody- 
namic limit, i.e it becomes independent on the length l x 
of the chain. Thus we have in the thermodynamic limit 



Si < l y c x . 



(18) 



To obtain an upper bound to S 2 we use the inequality 
F{x) < 1 - In 2 + In a;. This yields with eq. ([15]) 



2a / 2x 

S 2 <Ul-ln2) + -£lnk£ 
k=i \ j=i 



w kj aj(Ao ■ Di) 



To further evaluate the last term we make use of the 
convexity of the logarithm together with the arithmetic 
mean inequality 



fe=i \ i=i 



w k jaj(Ao ■ Di) 



I htJ^J^^a^Ao-Di) ] (20) 
3=1 k=i 



■In f£ ai (Ao-Di) 



3=1 



where we have used = 1 in the last step. 

We now have to evaluate the remaining logarithm. For 
this we make use of the fact that A and Q are regu- 
lar (i.e. strictly positive) Toeplitz matrices. Because 
of this, their elements can be obtained from_the non- 
negative spectral functions A (9) and q (9) [ll| Ak = 

i / A (6>) exp [-ik9] d9 , 


Qk = -i= J q (9) exp [— ik6\ d9. Since we have assumed 


above that also A — Q is strictly positive, the functions 
A (9) , q(9) are strictly positive and A (8) > q(9). In 
addition, we require also that (A (9) — q {9)) ±1 and q (9) 
have bounded derivatives of second order. As a conse- 
quence one finds (see [ll| page 221) 

which is a constant independent on l x . Thus the desired 
upper bound to the entropy for sufficiently large l x , l y is: 



Si = (ly - 1) E F (VMAo-Do] 



(17) 



fc=i 



Since A — Q was assumed to be strictly positive, the in- 
teraction V has only nonzero eigenvalues and thus corre- 
sponds to a gaped oscillator chain. As shown in the 



lx 

S < c\l y + c 2 l x + — lnl y 



(22) 



where c\ , c 2 are some constants independent of the size 
of the subsystem. 

A lower bound to the entropy can be found from the 
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inequality F (x) > In a;. This yields, with eg. (1151) . 



s ; 



In [a k (A ■ D 



(23) 



fe=i 



+|ln(g + ^ln (^^a i (Ao- J Di)) . 

fc=i y=i / 

Making use of Jensen's inequality for concave functions 
m \ t j a j) ^ J2j tj l n (aj) and J2k w kj = 1 we find 



S > 



(ly ~ 1) 



Jin (a k (Ao-Do)) 



(24) 



k=l 



j'=l 

To evaluate the sums over the logarithms we employ 
Szego's theorem [ll[ for determinants of a Tocplitz ma- 
trices T. The theorem states: for sufficiently large l x 

oo 

ln(dct {T)) ^ q Q l x + Y,k\q k \ 2 , 
fc=l 

for regular spectral function <?(#). Here % is Fourier 
coefficients of In q (8) . Since moreover 



£3 In (^-(Ao -.Di)) = ln( JJ ai (i4o-Di) 

= ln[dct(A )det(Di)], (25) 
we eventually find the lower bound 



5 > ai/ x + a 2 l y + -^ln(l y ). 



(26) 



Here a\,a 2 are constants independent of the size of the 
subsystem and we have ignored an unimportant constant 
term. 



By combining the two estimates (|22j) and (|26|) one finds 



c\ l x + C2 l y + ln(/ y ) > S > axl x + a 2 l y + -^\a(l y ). 

Since both sides of this inequality have the same func- 
tional form, S approaches for large l x ,l y the asymptotic 
value 



S 



: ln(Z y ), l x , l y > 1. 



(27) 



This is the main result of our paper, ft shows that the 
entropy area law is violated for a set of harmonic chains, 
which for themselves have a gaped spectrum and are non- 
critical but become gapless by a collective interaction be- 
tween the chains. Both upper and lower bound to the 
entropy attain the same logarithmic correction term to 
the area law. 

A physical system that can be approximated by the 
model studied here is an anisotropic ion crystal. In such 
a system the Coulomb-interaction in the direction of the 
small lattice constant can in first approximation be con- 
sidered as collective, while the one in an orthogonal di- 
rection is of finite range. 

In conclusion, we derived an exact asymptotic expres- 
sion for the entanglement entropy of a critical system of 
interacting oscillators in more than one dimension. We 
found that similar to one-dimensional systems @ the en- 
tanglement area law is violated by a logarithmic correc- 
tion proportional to the surface area in the critical di- 
rection. To our knowledge the system of collectively in- 
teracting harmonic strings considered here, which is ap- 
proximately realized e.g. in an anisotropic ion crystal, is 
the first nontrivial example of a critical two-dimensional 
system for which the correction to the area law can ex- 
plicitly be calculated. 
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